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DEVICE AND METHOD FOR ANALYZING A SAMPLE PLATE 
Background Information 

The present invention relates to a device for the analysis of 
a sample plate according to the type more closely specified in 
the generic part of Claim 1, as well as a method for analyzing 
5 the sample plate. 

It is of major importance, particularly in the field of 
material sciences, chemistry and pharmacy, to discover and 
develop optimized substances and materials with respect to 
respective cases of application. In this connection, sensor 

10 systems represent a special field of application, which are a 
key technology having a steadily growing number of 
applications, both in industry and in the private sector. For, 
sensors are used, for example, in technical process 
monitoring, in the area of environmental protection, in the 

15 field of medicine and in the motor vehicle field. A 

considerable quantity of development work is currently being 
put especially into the development of faster and highly 
sensitive sensors having a low cross-sensitivity. 

As a rule, up until now, the developments were limited to 
2 0 optimization or modification of known materials. However, 

there is the problem that, for certain areas of application of 
the sensor systems, there is a great requirement for new 
materials which cannot be sufficiently covered by conventional 
methods which are distinguished by the production of 
2 5 individual sensors and a subsequent sequential 
characterization . 

In particular, in the development of new types of sensitive 
materials or material combinations, it may be expedient to use 
methods from the field of combinatorial chemistry or so-called 



NY01 1069865 vl 



high throughput methods. In these methods, parallelized 
synthesis methods and screening methods are involved, by which 
new materials and material combinations may be discovered, or 
known synthesis methods for existing materials may be 
5 optimized in a broad field of parameters. 

A general representation of high throughput methods is known 
from US 5,925,355, and in this document, especially the 
application of combinatorial chemistry, known essentially from 
the field of pharmacy, to chemical and material science 
10 application areas being proposed. 

A device of the type named in the introduction, for analyzing 
a sample plate is known, for example, from DE 101 31 581 Al. 
This device includes a sample plate on which 64 material 
samples have been applied in matrix fashion, which, in each 
15 case, are connected to two electrodes that, in turn, are 
provided with contact locations to which a means for 
reversible and addressable contacting to the installation may 
be brought . 

Summary of the Invention 

20 The device, according to the present invention, for analyzing 
a sample plate having the features according to the generic 
part of Claim 1, and having a measuring head, that may be 
positioned in a housing carrier, which includes two measuring 
wires, per material sample for the electrical connection to 

25 the contact means, which lie against contact surfaces of the 
sample plate in a prestressed manner, and are connected to a 
measuring unit and an evaluation unit, has the advantage of 
being simple to handle, since, in the use of a normalized 
sample plate, the contact between the material samples and the 

3 0 measuring unit and evaluation unit is able to be produced by 
simply setting the measuring head into the housing carrier. 
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The device according to the present invention is particularly 
suitable for developing and discovering materials or material 
combinations that may be used as sensor materials, and may be 
characterized by their electrical properties. For example, the 
5 device according to the present invention may be used for 
developing an optimized sensor material of a gas sensor. 

Using the device according to the present invention, it is 
possible almost simultaneously to investigate a large number 
of potential sensor materials, disposed on the sample plate, 
10 under various test gases at different temperatures which may 
amount to up to 800°C. The investigation may be made by 
potent iometric, resistive, capacitative methods and also by 
complex impedance spectroscopy. 

The concept of spectroscopy should be understood, in this 
15 case, as involving frequency-dependent measurements, that is, 
the impedance of a sample is investigated at different 
measuring frequencies. For example, the individual material 
samples are tested in each case in a frequency range between 
10 Hz and 10 7 Hz at a measuring data density of 15 measuring 
20 points per decade. This means that, per material sample, 180 
measured data are ascertained at such a measuring data 
density. A far-reaching data reduction may be achieved, for 
instance, by the adaptation of a suitable circuit equivalent 
to the measured data . 

25 In order to hold the measuring wires which, in particular, lie 
against the contact surfaces of the sample plate via a fusion 
ball, with prestressing, the measuring wires may each be 
connected to an especially gilded spring contact, which will 
ensure a constant contact pressure of the respective measuring 

30 wire against the respective contact surface. 
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In order to be able to apply a certain test gas to the 
material samples, for the development of a gas sensor, the 
measuring head may be connected to a gas supply unit. 

In order to be able to expose the material samples to 
5 different test gas atmospheres or reference gas atmospheres, 
the gas supply unit, which is expediently connected to a data 
processing unit of the measuring and evaluation unit, includes 
a gas mixing device. Furthermore, the gas supply unit may 
include a water reservoir for moistening the test gas 
10 atmosphere or reference gas atmosphere. 

Furthermore, the measuring head may be designed in such a way 
that it includes a gas chamber, as an integrated component, 
that is situated above the material samples of the sample 
plate and is preferably formed by an essentially bell -shaped 
15 distributing device. The gas chamber is connected to the gas 
supply unit . 

In order to achieve a homogeneous distribution of the test gas 
or the reference gas in the gas chamber, in one advantageous 
specific embodiment of the device according to the present 

2 0 invention, a diffuser is disposed in the gas chamber. 

In the case of a plurality of material samples on the sample 
plate, in order to be able to measure the individual material 
samples in a simple manner, the measuring and evaluation unit 
advantageously includes two relay switch panels that are 
25 connected to the measuring wires, and, for example, in the 
case of 64 material samples on the sample plate, preferably 
each have a 3x64 matrix of relays appropriate for high 
frequency. In this case, the 64 material samples on the sample 
plate are able to be measured in a measuring cycle, at least 

3 0 three measuring variables being accessible, namely, for 

example, the impedance via an impedance analyzer and direct 
current resistances of the material samples via additional 
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appropriate measuring devices and their current /voltage 
characteristics curves . 

The measuring and evaluation unit is expediently equipped with 
a measuring and control software, which, on the one hand, 
5 controls the measuring progress and, on the other hand, passes 
along the measured data obtained to an appropriate data file 
or even to a relational databank, which may be read in by an 
evaluation software . 

The evaluation software preferably works in such a way that it 
10 includes a fit functionality for computing theoretical 

impedance spectra for the individual samples, the computation 
preferably being made based on a circuit equivalent that 
includes at least one virtual or real electronic component- A 
virtual component is, for instance, a constant phase element 
15 (CPE) . Thus, the fit functionality computes, based on the 

measured data for a circuit equivalent made up, for example, 
of a serial RC element, a theoretical impedance spectrum 
approximated as well as possible to a measured impedance 
spectrum, a variation of the capacitance and/or the resistance 
2 0 of the components of the RC element being carried out for the 
adaptation of the theoretical to the measured spectrum. Thus, 
in this case, a resistance value and a capacitance value for 
the RC element are assigned to the theoretical impedance 
spectrum that is adapted as well as possible. From the 

2 5 resistance value, for example, one may draw conclusions on the 

sensitivity of the respective material sample by comparing it 
to a reference value. 

And, in the case of a plurality of material samples, which are 
measured under various measuring conditions, in order to make 

3 0 accessible to a simple evaluation the output values obtained 

for the individual material samples representing, for 
instance, the sensitivities, the evaluation software 
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advantageously includes a datamining functionality. The 
datamining functionality determines the optimal material 
sample for the respective application case by a numerical 
path, for instance, by advantageously applying 
5 multidimensional object functions. 

Alternatively or in addition, the datamining functionality may 
also include a visualization functionality. In this case, a 
user is able to determine the optimal material sample for the 
respective application case, supported by a video screen. For 
10 instance, the visualization functionality works with a color 
spectrum, one certain color having a high sensitivity of the 
material sample assigned to it and another color having a 
lower sensitivity assigned to it. 

In order to be able also to discover sensor materials that are 
15 suitable for application cases in which, such as for instance 
in exhaust gas, where high temperatures prevail, the device 
according to the present invention preferably has a heating 
device into which the sample plate may preferably be immersed. 
In this case, the device especially has an high temperature 
20 reactor that is bordered by the heating device, and in which 
the sample plate that includes the material sample may have 
applied to it different test gases or reference gases. 

The present invention also has a method for analyzing a sample 
plate having the features of Claim 19 as subject matter. By 

25 applying this method, it is possible electrically to measure, 
under different conditions, a sample plate on which a large 
number of material samples, for instance, 64 material samples 
are disposed, and to select in a fully automatic fashion a 
material sample best suited to the application case. This 

30 takes place by an automatic selection of starting values for 
the components of the respective circuit equivalents and the 
subsequent error minimization computation. The starting values 
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are applied in the error minimization computation, and 
starting from the starting values, a theoretical impedance 
spectrum, adapted to the respectively measured impedance 
spectrum, for the respective material sample being calculated 
under the respective measuring conditions. By this procedure, 
a considerable data reduction is made possible, since, 
starting from the impedance spectra which have a plurality of 
measuring points, few fit values representing derived 
variables are ascertained that describe the individual 
components of the circuit equivalents. The fit values thus 
represent dimensionings of the components of the circuit 
equivalents, using which the measured impedance spectrum may 
be simulated as best as possible. 

In the determination of the circuit equivalents, for example, 
15 a serial connection of four RC elements may be selected, the 
higher RC elements, if necessary, being set to a starting 
value that has no influence on the computation of the 
theoretical impedance spectrum. 

The starting values of an RC element, required for the error 
20 minimization computation, are preferably computed from the 
maximally measured, imaginary impedance Z"_MAX and the 
corresponding measured frequency f_Z"_MAX, according to the 
following formulas : 

Rl START = -2Z"_MAX 

CI START = . 

In f _Z" _MAX R _START 

25 

If the circuit equivalent includes several RC elements, the 
starting values for the error minimization computation 
required for the higher RC elements are preferably ascertained 
from the difference spectra between the measured data and data 
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which are computed or simulated based on the starting values 
computed for the first RC element. 

In particular, the selection of so-called "good" starting 
values, which are close to the actual variables of the 
5 components of the circuit equivalent, decisively shortens the 
duration of the subsequent error minimization computation. 
"Bad" starting values, on the other hand, may lead to the 
error minimization computation, carried out based on the 
starting values, yielding meaningless values. 

10 In one preferred specific embodiment of the method according 
to the present invention, the circuit equivalent for the 
simulation of impedance spectra, based on fit values obtained 
in response to an error minimization computation, is made up 
of a serial separation of four RC elements. By the 

15 determination of starting values for the individual 

components, it is then determined, in the light of a threshold 
value, how many RC elements are taken into consideration in 
the simulation computations. The threshold value is preferably 
a value that may be preset by a user. The percentage ratio of 

20 the resistance of the first RC element to the resistance of 
the current RC element n is checked according to the formula 

RCn_START > Wert [%] • RC1_START. 

The variable "value" is the variable, changeable by the user, 
which, for instance is preset to 10%. If the argument is not 
25 satisfied, the starting values for the components of the 

respective RC element are set to values that have no influence 
on a simulation of the impedance spectra. These values are 
kept constant in the error minimization computation. 

Furthermore, in the method, the validation variable is 
3 0 determined which evaluates the agreement between the 
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calculated, theoretical impedance spectrum and the 
respectively assigned, measured impedance spectrum. 



The evaluation quantity represents the output value relevant 
to the respective analysis which, for example, in the 
5 determination of a sensor material for a gas sensor, reflects 
a sensitivity of the respective material sample. The 
sensitivity is a measure for the quality of a sensor. 

The sensitivity of a resistive gas sensor may be defined in 
various ways. If, for instance, one takes into consideration 
10 the direction of the change of the resistances during the 

application of gas, a sensitivity S may be expressed as the 
quotient of the resistance R_TEST under a test gas atmosphere 
and the resistance R_0 under reference conditions, as follows 

c RTEST . 

S = — - — ; for oxidizing gases 

S = + — ^ -® — • f or reducing gases 
R _ TEST B 

15 The sign of sensitivity S gives information with regard to the 
resistance change . 

Alternatively, the sensitivity may be described as a change in 
the resistances, namely, according to the formulas 

0 . R TEST-R 0 , 
" ~R TEST' ' f ° X oxldlzmg gases 

. A R 0-R TEST ^ 

S _ A = H — = — — — = ; for reducing gases. 

20 

These formulas yield sensitivities S_A between -1 and 1, that 
is, normalized sensitivities. 

It is also possible mutually to convert the sensitivities of 
the two definitions into each other. The sensitivity expressed 
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by the value S is especially meaningful in the case of a large 
change of the resistance of a material sample that results 
from the application of a test gas. The sensitivity expressed 
by S_A is particularly meaningful if the resistance of a 
5 material sample changes only a little as a result of the 

application of a test gas. However, the sensitivity expressed 
by S_A has the advantage of a large tolerance with regard to 
measuring accuracies, a precise illumination of a region of 
small resistance changes, and therewith a better valuation of 
10 cross-sensitivities, as well as the possibility of an 
automated visualization and data processing. 

In one particularly simple specific embodiment of the method 
according to the present invention, the circuit equivalent is 
made up of a virtual arrangement of real electronic components 
15 such as capacitors and resistors. In its simplest embodiment, 
the circuit equivalent is composed of one resistor. 

By the error minimization computation, in which an adaptation 
of the variables of the components to the measured data takes 
place, a simulated impedance spectrum is determined for the 

20 circuit equivalent for the respective material sample. From 
the adapted variables of the components, one is able to draw 
conclusions on the electrical properties of the processes in 
the material sample described by the circuit equivalent. If 
the electrical behavior of the material sample is determined 

25 by several processes of different relaxation times, it is 

necessary to draw upon complicated circuit equivalents which, 
for instance, are made up of several RC elements connected in 
series. If the processes react differently to a variation of 
the measuring conditions, one may undertake an assignment of 

30 individual processes to components or groups of components so 
that the individual processes may be analyzed separately. 
Because of the error minimization computation, the data volume 
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for the description of the impedance measurement is reduced to 
the arrangement and the magnitudes of the components. 

In the case of a sample plate having, for example, 64 material 
samples and a measurement under eleven different gas 
atmospheres at four different temperatures in each case, the 
quantity of measured data is approximately 2816 impedance 
spectra. In the method according to the present invention, the 
individual impedance spectra are simulated and in each case 
reduced to the derived measured variables that are reproduced 
by the variables of the components of the circuit equivalent. 

The computation of a theoretical impedance spectrum of a 
circuit equivalent is made by one preferred specific 
embodiment of the method according to the present invention in 
such a way that a complex admittance Y* and a phase shift q> of 
the individual components of the RC elements are determined at 
a given angular frequency co (211 measuring frequency) according 
to the following formulas: 



0=arctan (coRC) 
Y' = cos0 
Y"= sin0 

The respective impedances may be assessed by a transformation: 




r 



T 



r 2 +y" 



where Z' is the real part of the impedance and Z" is the 
imaginary part of the impedance . 
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By a serial connection of the RC elements, the impedances of 
the individual RC elements may be summed directly. A 
calculation of the impedances for frequencies that correspond 
to the measured frequencies yields a data set corresponding to 
5 the measured data, so that an error estimation between the 
measured data set and the theoretically determined data set 
is possible. 

The error minimization computation carried out in the method 
according to the present invention is carried out in one 
advantageous embodiment by variation of the variables of the 
individual components by 1%. An error may be determined by 
analysis of the differences between the theoretically 
calculated spectra and the measured ones. If, after a 
variation, an error decreases, a renewed variation is carried 
out for the same component of the circuit equivalent, and, 
based on this variation, a theoretical impedance spectrum is 
computed, using which a new error calculation is carried out. 
If the error does not decrease, the variation that has taken 
place is reversed, and the variation of the components having 
the opposite sign is varied, or another component is varied. 

The error computation is preferably carried out in such a way 
that the function for determining the error takes into account 
the starting value of the resistance of the first RC element. 
If the value is greater than a nominal measured resistance of 

2 5 the impedance analyzer used in the measurement, for instance, 

greater than 3- 10 7 Q, the error is determined only from the 
imaginary part of the calculated impedance. In order to be 
able to weight the high-frequency range of the ascertained 
spectrum, the ascertained deviations are preferably multiplied 

3 0 by the logarithm of the measured frequency. This weighting 

makes it possible to suppress physically meaningless or even 
erroneous measured data in the low- frequency range. 

12 
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For the calculation of the validation variable, that is, for 
the validation of the quality of the theoretically calculated 
impedance spectrum for the respective material sample, 
preferably a corridor is determined about the theoretical • 
5 impedance spectrum which contains, for example, 90% of the 

measured data. In order to minimize the time required for the 
computation of the validation variable, a successive 
approximation algorithm may be used for this, the boundary 
values of the algorithm 0 Q and of the double values of the 
10 real parts of the summed impedance may be at the smallest 
measured frequency . 

The method according to the present invention is preferably 
used for the analysis of material samples under various test 
gas atmospheres and especially at various temperatures. In 
15 this case, impedance spectra are measured under various 

measuring conditions for all material samples situated on the 
sample plate. 

In the method according to the present invention, the 
respectively measured spectra are then each simulated in a 
2 0 theoretical manner by an error minimization computation based 
on a circuit equivalent. This results in a large quantity of 
evaluation magnitudes which represent the target variables in 
the method. 

In the method according to the present invention, the target 
25 variables or the evaluation variables are preferably written 
into a databank and evaluated using a datamining 
functionality. 

As far as the databank is concerned, which is used for 
receiving and making available data sets, expediently a 
30 relational databank is involved in which the data are filed in 
the form of tables, ordered according to thematic areas. 
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Relationships between the individual data sets in the tables 
are produced by so-called identification keys. 

The databank may, for example, include additional properties 
of the material samples, such as the synthesis conditions of 
5 their starting materials, their sampling history and the like. 
These properties are linked with one another via table 
relationships . 

The datamining may be performed using a possibly 
multidimensional target function and/or using a visual 
datamining functionality. Datamining, carried out using a 
target function, is a numerical method which is based on the 
individual evaluation variables which are stored in the 
databank, for instance, classified according to measuring 
temperature and test gas. In this context, it is first 
specified which properties are wanted of the material being 
sought. For instance, in the search for a sensor material for 
a gas sensor it may be stated with respect to which gas the 
sensor material is to be sensitive, and which cross- 
sensitivities could have an interfering effect. Thus, from 
this there develops a demand profile for a fingerprint of the 
sensitivities . 

The visual datamining functionality advantageously works in 
such a way that the evaluation variables of the material 
samples of the sample plate are represented broken down, for 
25 instance, according to test gases and temperatures. 

In the use of the device according to the present invention 
and the method according to the present invention, one has a 
preferably fully automatic high- throughput impedance system 
using which new materials having a high sample throughput and 
3 0 a low time expenditure and cost expenditure are able to be 

developed. Thus, for example, by using the system on two days, 
64 different material samples may be investigated and also 
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evaluated at four different temperatures and under eleven 
different test gas atmospheres, with regard to their sensoric 
properties . 



In particular, the system may be used for general material 
5 development and especially for the development of sensors in 
the motor vehicle field and in the safety technology field. 

The present invention also has data processing equipment 
having a data processing program as subject matter, for 
implementing the method according to the present invention. 

10 Further advantages and advantageous refinements of the subject 
matter in accordance with the present invention result from 
the description, the drawing and the patent claims. 

Brief Description of the Drawings 

An exemplary embodiment of the device according to the present 
15 invention is shown in a simplified version in the drawing and 
is elucidated in more detail in the following description, in 
connection with a method according to the present invention. . 
The figures show: 

Figure 1 a test setup of a device according to the present 
20 invention; 

Figure 2 a measuring appliance of the device according to 
Figure 1 ; 

Figure 3 a measuring head of the measuring device according 
to Figure 2 ; 

25 Figure 4 a gas distribution device of the measuring head as 
in Figure 3 ; 

Figure 5 a measuring sequence in the light of a flow diagram; 
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Figure 



a diagram representing a measuring sequence over 
time for material samples at various temperatures 
and under various gas atmospheres; 



Figure 7 a diagram in which an influent behavior is shown of 
5 four different surface doped material samples; 

Figure 8 a diagram in which a response time performance is 
shown of three different material samples during a 
test gas pulse; 

Figure 9 a circuit equivalent; 

10 Figure 10 a flow diagram of an error minimization computation; 

Figure 11 an exemplary curve of the sensitivity of a material 
sample; 

Figure 12 an exemplary fingerprint of sensitivities of a 
material sample with regard to various gas 
15 atmospheres, and 

Figure 13 a greatly schematized representation of a visual 
datamining f unct ional i ty . 

Detailed Description of the Exemplary Embodiment 

In Figures 1 through 4, a test setup is shown of a device 10, 
2 0 for analyzing a sample plate 12, on which 64 material samples 
13 are situated. Device 10 includes a high temperature reactor 
14, a gas supply unit 16 and a measuring and evaluation unit 
18. 

High temperature reactor 14, which is shown in more detail 
25 especially in Figures 2 through 4, includes a frame or housing 
20, in which an adjustable-height heating block 22 is 
situated, which is supported on a guide rod 23 and a threaded 
rod 24 that is provided with a crank 25. 
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Heating block 22 , which is made up, for instance, of four 
heating plates that border on a heating chamber, having a 
heating power of 1100 W each, is adjustable in height in such 
a way that plate carrier 28, and thus sample plate 12 may be 
5 immersed into the heating chamber. 

Furthermore, high temperature reactor 14 includes a measuring 
head 26 which is set into measuring head carrier 27 that is 
connected to frame 20, and is used for contacting the 64 
material samples 13 situated on sample plate 12 to measuring 
10 and evaluation unit 18. Sample plate 12 is situated on a plate 
carrier 28 that is connected to measuring head carrier 27. 

Measuring head 26, which is particularly shown in Figure 3, 
includes a ground plate 29, which is formed, for example, from 
a workable glass ceramic and is used as a mounting support of 

15 128 measuring wires 30A, 3 OB, that are each made of platinum 
and are sheathed by an aluminum oxide tube. At their lower 
ends, measuring wires 3 OA, 3 OB each have a fusion ball 31A, 
3 IB which are used for contacting each measuring wire to a 
contact surface of sample plate 12 . Sample plate 12 has two 

20 contact surfaces per material sample, i.e., in the present 

case altogether 128 contact surfaces, which each collaborate 
with one of measuring wires 3 OA, 3 0B. Measuring wires 3 OA, 3 0B 
are provided at their upper ends, in the region of ground 
plate 29, each with a gilded spring contact 32A, 32B, which, 

25 when measuring head 26 is in a fixed position, ensures a 

constant contact pressure of measuring wires 3 OA, 3 0B and of 
fusion ball 31A, 31B against respectively assigned* contact 
surfaces of sample plate 12. Furthermore, the sheathings of 
measuring wires 30A, 30B, that are made of platinum, are fixed 

3 0 to supporting plates 33 and 34, which are fastened to a rod 
assembly 35, and are aligned parallel to ground plate 29. 
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Spring contacts 32A, 32B are connected via lines 36A, 36B to 
SMB sockets 37A, 37B that are integrated into measuring head 
26, and, in turn, are connected via shielded SMB lines 38A, 
38B to measuring and evaluation unit 18 of device 10. 
5 Altogether, measuring head 26 has 128 SMB sockets 37A, 37B, of 
which each is connected to a measuring wire 3 OA, 3 OB and to 
which, in each case, one SMB line 38A, 38B is connected, that 
leads to measuring an evaluation unit 18. However, for 
clarity's sake, in Figures 1 to 3 , in each case, only two SMB 
10 sockets 3 7A, 3 7B and two SMB lines 3 8A, 3 8B are shown, as well 
as measuring wires 30A, 3 OB that are allocated in each case. 

Furthermore, measuring head 26 has a gas distribution device 
39, made, for example, of quartz glass and essentially bell- 
shaped, which is shown in detail in Figure 4 and which is 
15 connected via a gas supply line 40, made of high-grade steel, 
to gas supply unit 18. Bell-shaped gas distribution device 39 
borders on a gas chamber that is situated above the 64 
material samples 13 that are disposed on sample plate 12. 

Material samples 13 are formed, for instance, of tin oxide 
20 Sn0 2 , and have different dopings which are formed, for example, 
of lanthanides. Material samples 13 are distributed in a 
matrix-like manner, in eight rows and eight columns, on sample 
plate 12 . 

Gas distribution device 39, that is shown in detail in Figure 
25 4, also has in its gas chamber 41 a diffuser insert 42, formed 
of a quartz ball and has a plurality of bore holes 43 that 
each have a diameter of about 1 mm. 

In order to manage a uniform gas outflow from gas chamber 41, 
bell -shaped gas distribution device 39 has at its edges 
30 spacers 43, which establish a gap of 0.8 mm in width between 
gas distribution device 39 and sample plate 12. 
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In order for gas chamber 41 , and therewith material samples 
13 , to have different test gas atmospheres applied to them, 
gas supply unit 16 has two gas bottle cabinets, not shown here 
in detail, each having four gas bottles, which are 
5 respectively connected to a gas flow controller 44A, 44B, 44C, 
44D, 44E, 44F, 44G and 44H, one gas bottle containing moist 
synthetic air, the second gas bottle containing hydrogen, the 
third gas bottle containing methane, the fourth gas bottle 
containing synthetic air, the fifth gas bottle containing 

10 nitrogen dioxide, the sixth gas bottle containing nitrogen 
monoxide, the seventh gas bottle containing propane and the 
eighth gas bottle containing carbon monoxide. The capacities 
of flow controllers 44A, 44B, 44C, 44D, 44E and 44F are 
respectively between 0 seem and 100 seem. The capacities of 

15 flow controllers 44G and 44H are respectively between 0 seem 
and 10 seem. By appropriate control of the volume flow of the 
various gases, using the group of uniform devices, coming from 
the eight gas flow controllers 44A, 44B, 44C, 44D, 44E, 44F, 
44G and 44H test gases of different compositions may be fed 

20 into a collecting line 45 that is connected to gas supply line 
40. 

In order to set a relative humidity of the respective test 
gases, the test gas may have admixed to it a moist carrier 
gas, consisting, for example, of synthetic air. The moisture 
25 in the carrier gas is set by bubbling it through a water 

reservoir 46. Measuring the humidified carrier gas using a 
moisture sensor, not shown here in greater detail, yields, for 
example, a relative humidity of approximately 90% at room 
temperature . 

3 0 Measuring and evaluation unit 18 includes two relay switch 

panels 50 and 51, to which respectively 64 lines 38A, 38B are 
connected that lead to measuring head 2 6 of high temperature 
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reactor 14. Relay switch panels 50 and 51 each form a 3x64 
matrix of relays suitable for high-frequency. 

Relay switch panels 50 and 51 are connected to a measuring and 
evaluation computer 53 via a digital control line 52. 
5 Moreover, relay switch panels 50 and 51 are connected to an 
impedance analyzer 64 and a so-called source meter 55, via 
measuring lines 54. These two measuring unite are also 
connected to measuring and evaluation computer 53 via digital 
control line 52. The addressing of impedance analyzer 54 and 
10 source meter 55 takes place via the two relay switch panels 50 
and 5 1 . 

Measuring and evaluation computer 53 is also connected via an 
additional digital control line 56 to a D/A-A/D converter 57, 
which is connected via an analog control line 58, on the one 
15 hand, to heating device 22 of high temperature reactor 14, 

and, on the other hand, to gas flow controllers 44A, 44B, 44C, 
44D, 44E, 44F, 44G and 44H of the gas supply unit or gas 
mixing battery 16. 

On measuring and evaluation computer 53, a modular measuring 
20 and control software is filed, which makes possible, via a 

script control, a complete automation of measurements carried 
out using device 10. Because of the modular design of the 
measuring software and control software, a broadening of the 
measurement system is possible without a problem. 

25 Using source meter 55, direct current resistances, U/I 

characteristics lines or voltages of individual material 
samples 13 may be measured. These values, as well as the 
measured values ascertained using impedance analyzer 54, may 
be passed on via digital control line 52 to measuring and 

30 evaluation computer 53, which includes a databank for the 
measured data. 
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In order to measure material samples 13, that are disposed on 
sample plate 12, in a high throughput mode or according to a 
high throughput method, via a script data file is used, a list 
of tasks that is made up of index key words and parameters 
5 being passed over, via a script data file, to the software 
using which the script data file is processed and which 
controls all the functions of the system. The script data 
files ensure a continual checking of control parameters, so 
that measurements under faulty measurement conditions are 
10 excluded in that a further processing of the script data file 
is suspended. There is no limit with regard to the length of 
the script data file. 

A measuring sequence able to be carried out using the device 
according to Figure 1 is shown in Figure 5 in the light of a 
15 flow chart. For the analysis of resistive gas sensor 

properties, material samples 13 of sample plate 12 are first 
electrically characterized under a reference gas atmosphere 
which, for example, is formed of synthetic air of a relative 
humidity such as 45%, and under various test gas atmospheres. 

2 0 By a modulation of the measuring temperature, data may be 

gathered on the influence of the respective operational 
temperature on the material samples that are suitable as 
sensor material, as well as on the activating energy of 
conductivity processes. 

25 In order to ensure comparability of measured data of several 
sample plates, the script data file used is designed as 
standard script which controls the complete high throughput 
screening, whose chronological sequence is described in the 
flow chart shown in Figure 5. In this instance, in a first 

3 0 method step Ml, a measuring temperature T_MESS of sample plate 

12 is set. Measurements may be carried out, for example, in a 
temperature range between 400° C and 250° C, in steps of 50°. 
During the cooling-off phases, sample plate 12 has applied to 
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it a reference gas of synthetic air having a relative humidity 
of 45%. 

When there is a change in the sample temperature, a stable 
base resistance of the material samples first sets in after a 
5 certain time, so that conditioning of the material samples in 
a step M2 is required. Temperature changes in the material 
samples lead to a metastable state of intrinsic defects, which 
may represent oxygen vacancies, and whose thermodynamic 
equilibrium setting may require a finite time. Figure 7 shows 

10 in an exemplary manner the conditioning behavior or influent 
behavior for four different surface-doped ln 2 0 3 samples upon 
reaching a target temperature of 300° C. Resistances R of the 
samples approach a boundary value, asymptotically with time, 
which represents the so-called base resistance or reference 

15 resistance. This boundary value is reached, as extrapolated, 
after about 90 minutes. To ensure a constant reference 
resistance, it is of advantage to select a conditioning time 
of 120 minutes. 

Subsequently, the test gas required for the measurement is put 
20 together in a method step M3 , and is introduced into gas 

chamber 41 above sample plate 12. The first test gas includes, 
for instance, hydrogen at a concentration of 2 5 ppm. Synthetic 
air acts as the carrier gas. The gas flow is set to 100 seem 
using moist synthetic air. 

2 5 In order to ensure that the individual material samples have 

reached their base resistance independent of their relative 
position on the sample plate, there subsequently takes place, 
in a step M4, an advance phase for the test gas that was 
analyzed before. Figure 8 shows the resistance pattern of 

3 0 three material samples lying on a diagonal of the sample 

plate, during a test gas pulse lasting 40 minutes, the pulse 
gas containing propane at a concentration of 50 ppm. The base 
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material of the material samples is tin oxide Sn0 2 in the 
example shown in Figure 8, because of its great sensitivity to 
hydrocarbons. Independently of the relative position on the 
sample plate, the resistance of the material samples 
decreases, because of the test gas pulse, to a constant value, 
within approximately 6 minutes. Thus, the response time 
amounts to 6 minutes in each case, the response behavior being 
essentially independent of the position of the respective 
material sample on the sample plate. After terminating the 
test gas pulse, the resistances tend asymptotically to the 
base resistances within approximately 10 minutes. Accordingly, 
the response behavior is essentially independent of the 
position of the material sample on the sample plate. The 
advance time, selected for the script data file, of each test 
gas or reference gas amounts to about 15 minutes before the 
corresponding measurement . 

In a subsequent method step M5 , there takes place the 
measurement of impedance spectra for the 64 material samples 
disposed on the sample plate. In this context, in the script 
2 0 data file, for the measurements of the impedance spectra, the 



following are established as parameters: 

amplitude of the measuring voltage [V] 0.1 

starting frequency [Hz] 10 

ending frequency [Hz] 10 7 

25 measuring points per frequency decade 15 

bias [V] 0 

mode [HS : highs peed; NO: normal; AV: average] HS 



The measuring data obtained are ascertained using impedance 
analyzer 64 and passed on to measuring and evaluation computer 
30 53 or the databank filed on it. 

Thereupon the measurement may be carried out under another 
test gas atmosphere, in this instance, once again the program 
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runs through steps M4 and M5 . For example, as further test 
gases, carbon monoxide at a concentration of 50 ppm, nitrogen 
monoxide at a concentration of 5 ppm, nitrogen dioxide at a 
concentration of 5 ppm or propane at a concentration of 25 ppm 
5 may be used. 

The entire screening may then be carried out by jumping back 
to step Ml at another measuring temperature . 

Figure 6 visualizes a measuring sequence determined by a 
standard script, in this diagram, measurements under reference 
conditions and during conditioning phases not being shown, 
for the sake of clarity. Measurements under reference 
conditions, i.e. under a reference gas atmosphere, are always 
carried out before the initiation of a new test gas H 2 , CO, NO, 
N0 2/ or propane at the appropriate temperature T. There is 
always respectively a so-called test gas advance phase X and a 
measuring phase Y. Furthermore, the respective concentrations 
C of the test gases are shown in Figure 6 . 

The measuring and evaluation software may have a functionality 
on which the respectively measured impedance spectrum is 
20 graphically shown on a monitor. The measured data, which are 
preferably stored in ASCII format in the databank or in the 
index assigned to the sample plates on the measuring and 
evaluation computer, are consequently directly accessible to 
visual control . 

2 5 Furthermore, the measuring and evaluation software may be 
furnished with a functionality which shows the impedance 
spectra of the sample materials of a sample plate as an image 
matrix independently of position. In this context, raw data as 
well as derived data may be shown. The functionality may also 

30 include an evaluation window to which the data of a certain 

material sample may be transmitted. The functionality may also 
be combined with additional image -producing measuring or 
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evaluation systems. Image data from these systems may then 
also be read into the image matrix independent of position, in 
order thus to be able to obtain further data on the samples . 

After measuring the impedance spectra for the individual 
5 material samples, a theoretical impedance spectrum based on 
circuit equivalents is computed for each material sample. In 
the present case, the circuit equivalent is made up of a 
serial connection of four RC elements of the kind shown in 
Figure 9. The theoretically calculated impedance spectra, in 
10 this context, are each adapted as well as possible to the 
correspondingly measured impedance spectrum, namely, by 
varying the magnitudes or dimensionings of the individual 
components of the circuit equivalents. 

This is clarified in the following with reference to the error 
15 minimization computation in the light of the flow chart shown 
in Figure 10. 

In order to carry out the error minimization computation, it 
is necessary to determine starting values for the individual 
components of the circuit equivalents. At the same time, it is 

20 necessary to determine the number of RC elements to be taken 
into consideration in the error minimization computation. For 
these purposes, in a step SI, the maximum imaginary impedance 
Z U _MAX is ascertained from the measured data for the 
respective material sample. In the selection, it is also 

2 5 judged whether, in the case of the maximum measured imaginary 
impedance, a fault measurement is involved, which is done by 
an investigation of the imaginary impedances measured at 
adjacent measuring points, i.e. by an investigation of the 
local maximum. 

30 Subsequently, in a step S2 , starting values R1_START, for the 
resistance, and C1_START, for the capacity of the first RC 
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element, are calculated based on Z"_MAX and the corresponding 
measuring frequency, according to formulas 

Rl _ STAR T = -2 • Z" _ MAX 

CI START = 

2/r • / _Z" _MAX • R _START 

5 where f_Z xx _MAX is the measuring frequency for the maximally 
measured, imaginary impedance. 

Subsequently, in a loop S3, starting values for the 
resistances and capacitances of the additional three RC 
elements of the circuit equivalents are ascertained. 

10 To do this, first of all, in a step S4 , based on starting 

values R1_START and C1_START, a theoretical impedance spectrum 
is calculated, and, in a step S5, a difference spectrum 
between the theoretical impedance spectrum and the measured 
impedance spectrum is calculated, and from this difference 

15 spectrum, in turn, the maximum of the imaginary impedance 

Z"_MAX is ascertained. Starting from this maximum, in a step 
S6, the starting values Rn_START and Cn_START (n = 2 to 4) are 
then calculated for the respective RC element being examined, 
according to the formulas explained in connection with step 

20 S2 . 

Thereafter, in a step S7, a threshold value examination is 
carried out for calculated resistance starting value Rn_START, 
using which, it is determined whether the respective RC 
element should be taken into consideration in the subsequent 
25 simulation calculations. The threshold value is a value that 
may be changed by the user. In the comparison, the percentage 
ratio of the resistance of the first RC element to the 
resistance of the current component is checked. 
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If the calculated starting value Rn_START is greater than the 
threshold value, the RC element is taken into consideration, 
the system goes back to step S3, and a calculation of the 
starting values Rn_START and Cn_START is undertaken for the 
5 next RC element. If the argument is not satisfied, the 

magnitudes of the components of the higher RC elements are set 
to values that have no influence on a simulation calculation 
for an impedance spectrum. In subsequent adaptation steps, 
these values remain constant, and they amount, for example, to 
10 Rn_START = 1 and Cn_START = 10" 15 . This determination takes 
place in a step S8 . 

Thus, the threshold value query determines the number of the 
RC elements that are to be taken into consideration in the 
following steps. If the threshold value is undershot, the 
15 system transfers directly to the error minimization 

computation based on the number of RC elements determined, for 
instance, based on an RC element, in this case, a loop S9 
being executed over all m RC elements taken into 
consideration . 

2 0 Within loop S9, an additional loop S10, run 1 to 3 is 

executed, in this case an empirical number of run-throughs 
being involved, which are carried out to increase the accuracy 
of the error minimization computation. 

Within loop S10, first of all, in a step Sll, resistance value 
2 5 R1_JSTART of the first RC element is drawn upon as variable 

VAR. The variation of the individual components in the error 

minimization computation amounts to 1%, which is expressed by 

a value FAK = 0.01 in step Sll. The error is preset to 10 99 . 

The variation of the respective component takes place in a 
30 step S12, according to formula VAR = V-FAK. A calculation of a 

theoretical impedance spectrum then takes place, based on the 

varied value, in a step S13 . 
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After that, an error is calculated which is based on a 
comparison between the impedance spectrum calculated for the 
correspondingly adapted component variables and the measured 
impedance spectrum. A function for determining the error is 
5 dependent on the starting value of the resistance of the first 
RC element. In order to calculate the error, in a step S14 it 
is checked whether the respective starting value is greater 
than the nominal measuring resistance of the impedance 
analyzer used, for example, greater than 3 x 10 7 Q. If this is 
10 so, the error is determined in a step S15, only by examining 

the imaginary parts of the impedance, according to the formula 

error =| Z" _ fit - Z" _ mess | •log/ 

for the weighting of the high-frequency range of the spectrum 
a multiplication being made by the logarithm of the measuring 
15 frequency, whereby physically meaningless or erroneous 

measured data in the low- frequency range are suppressed. 
Otherwise, in a step S16 an error examination is made, based 
on both the real parts and the imaginary parts of the 
impedances according to the formula 

20 error =\ Z" _fit-Z" _mess | + | Z' fit - Z' mess | . 

In a step 17, the system checks whether the error has become 
smaller because of the variation of the component in step S12 . 
If this is the case, the respective component, for instance, 
the resistance of the respective RC element, is varied once 

25 more by going back to step S12 . If the error has not gone 
down, the variation that has taken place in step S12 is 
reversed in a step S18, and, in a step S19 it is established 
whether the sign of the variation will be changed in a step 
S2 0, or whether, in a step S21 and a subsequent step 22, the 

3 0 next component, in this case the capacitance C, will be 

selected as the variable component. In this case, then, there 
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is once more a return to step S12 , the magnitude of the 
component is then again varied as long as the error is 
decreasing . 

The variation of the resistance and the variation of the 
5 capacitance are repeated by return to step S10 twice for each 
RC element. Thereafter, as a function of the number of the RC 
elements, that are to be taken into consideration, there is an 
n-fold return to step S9. 

If no further minimization of the error is possible by 
10 variation of the resistance and the capacitance of the RC 
element, the values for the resistance and the capacitance 
yielding the smallest error are output as evaluation variables 
of the error minimization computation. 

At the close of the error minimization computation taking 
place in steps S9 to S22, in a step 23, a validation function 
is ascertained for the theoretically determined impedance 
spectrum, by which the quality of the theoretically calculated 
spectrum may be estimated. Upon determination of the 
validation function, for the estimation of a quality of the 
error minimization computation a corridor is ascertained, 
around the calculated impedance spectrum, that includes 90% of 
the measured data. The error thus ascertained lies between 0 
and 1, an error of 0 expressing an ideal agreement of the 
measurement and the simulation calculation, and an error of 1 
expressing no agreement of the measurement and the simulation 
calculation. 

When using the method for establishing a sensor material of a 
gas sensor, sensitivities S_A are calculated, based on 
resistances of the circuit equivalents obtained by the error 
30 minimization computation. Sensitivities S_A thus obtained are 
normalized, and lie between -1 and +1, as may be seen in the 
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diagram in Figure 11. Figure 11 shows sensitivities for an 
assumed reference resistance of 100 Q, at varying test 
resistance R_test . Sensitivities S_A are stored in the 
databank. 

5 Subsequently, based on the ascertained sensitivities S_A, a 
numerical datamining is carried out, so as to be able to 
select an optimal gas sensor for a special application case. 
The datamining is carried out especially in a numerical way, 
the individual sensitivities S_A being accessed which are 

10 stored in the databank, classified by measuring temperature 
and test gas. It is first specified which properties are 
demanded of the desired sensor material. In the simplest case, 
it is stated for which test gas the sensor is to be used, and 
which cross-sensitivities will interfere. This yields a 

15 requirement profile on the fingerprint of the sensitivities, 
as shown in Figure 12 . The requirement profile is converted 
into so-called or requirements, whereby data sets 

having the desired properties are characterized. Additional 
conditions, such as the type of base material of the sensor 

2 0 material or its surface doping may be taken up into the target 
function. In the example in Figure 12, the requirement profile 
turns out so that all test gases A, B, C, E, except for test 
gas D, have a sensitivity less than 0.2 and greater than -0.2, 
and test gas D has a sensitivity greater than 0.9. In each 

25 case X states the reference. Using an evaluation 

functionality, the requirement profile is converted to an SQL 
filter scan and passed on directly to the databank. The 
results of the SQL instruction may then be shown in tabular 
form on a screen. 

30 Alternatively, or in supplementation, the sensitivities 

ascertained for the various material samples at the various 
measuring conditions may be evaluated using a visual 
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datamining functionality. In this functionality, which is 
shown in Figure 13, sensitivities of so-called library plates 
are shown classified as to test gases and temperatures on a 
screen. In the exemplary representation as in Figure 13, 
5 library plates 101 to 112 are shown for four temperatures A, 
B, C and D, and for three different test gases I, II and III. 
Each library plate 101 to 112 is assigned to one temperature 
A, B, C or D, and to one test gas I, II or III. Sensitivities 
S_A of the individual materials of library plates 101 to 112 

10 are each shown according to their position on the sample plate 
as circles 120 and in an off -shade. In Figure 13, for the sake 
of clarity, only four of the 64 material samples are shown for 
each library plate 101 to 112. Positive sensitivities S_A, for 
example, are shown in colors that go from black to yellow via 

15 red shades, whereas negative sensitivities are shown, for 

example, in colors that go from black to turquoise via blue 
tones. The diameter of individual circles 12 0 is determined by 
the validation of the error minimization computation, in the 
case of a large error, a small circle being shown, and in the 

2 0 case of a small error, a comparatively large circle being 

shown. Using this scheme makes an intuitive judgment of large 
data quantities possible for the user. 

The data sets required for the functionality shown in Figure 
13 are extracted directly from the databank, so that the 
25 measurements to be shown may be selected via filter functions. 
The number of the data sets indicated is only limited by the 
working memory of the evaluation computer used. 

Furthermore, it is possible in the visual datamining 
functionality that additional data, concerning the individual 

3 0 material samples, may be indicated by the screen- supported 

selection of a certain circle. 
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